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Using 116.1 fb 1 of data collected by the BABAR detector, we present an analysis of H° production 
in B decays and from the cc continuum, with the H° decaying into Q,~ K + and E~tt + final states. We 
measure the ratio of branching fractions B{E° Q,~ K + )/B(E° c -> E~tt + ) to be 0.294±0.018±0.016, 
where the first uncertainty is statistical and the second is systematic. The S° momentum spectrum 
is measured on and 40 MeV below the T(4S) resonance. From these spectra the branching fraction 
product B(B -» E°X) x B(E° -» H~tt + ) is measured to be (2.11 ± 0.19 ± 0.25) x 10" 4 , and the 
cross-section product a(e + e~ —* E^X) x B(E° C — > E~tv + ) from the continuum is measured to be 
(388 ± 39 ± 41) fb at a center-of-mass energy of 10.58 GeV. 

PACS numbers: 13.25.Hw,13.30.Eg,14.20.Lq 



In this Letter we present a study of the 2°(csd) [l| 
charmed baryon through two decay modes: S° — ► fl~ K + 
and S° — * S" ir + the former of which is expected to 
proceed almost entirely via internal W-exchange. We 
determine the ratio of branching fractions of these decay 
modes, which has been measured previously to be 0.50 ± 
0.21 ± 0.05 0, H It was predicted to be 0.32 with a 
quark model calculation in which no spin information is 
exchanged between quarks other than through a single 
W boson 

We also study S° production by measuring the spec- 
trum of the 2° momentum in the e + e~ center-of-mass 
frame (p*). A number of theoretical predictions for S c 
production in B decays have been made [1,0,H,@1- There 
are several possible production mechanisms, principally 
b — > ccs weak decays, b — > cud weak decays in which an 
ss pair is produced during fragmentation, and Cabibbo- 
suppressed b — > cus weak decays. At this point there 
is insufficient experimental evidence to determine which 
of these is the dominant mechanism, and no clear theo- 
retical consensus. Insight into the contributing processes 
can be gained by studying the shape of the p* spectrum. 
Evidence for S c production in B decays was presented 
previously by the CLEO collaboration, with a statistical 
significance of ~ 3er in the Sj? — * S~7r + decay mode and 
~ 4er in the S+ — ► S~7r + 7r + decay mode lOj. 

The data for this analysis were collected with the 
BABAR detector at the SLAC PEP-II asymmetric energy 
e + e~ collider; the detector is described in detail else- 
where ^lj- A total integrated luminosity of 116.1 fb -1 
is used, of which 105.4 fb -1 was collected at the T(45) 
resonance Q (corresponding to 116 million BB pairs) 
and 10.7 fb -1 was collected at a center-of-mass energy 
of 10.54 GeV, which is below the BB production thresh- 
old. These are referred to as the on-resonance and off- 
resonance data samples, respectively. 

The reconstruction of S° candidates takes place as fol- 
lows. A A candidate is reconstructed by identifying a 
proton and combining it with an oppositely charged track 
interpreted as a tt~ , fitting the tracks to a common ver- 
tex. The A candidate is then combined with a negatively 
charged track interpreted as a n~ (K~) to form a S" 
(f2~) candidate. For each intermediate hyperon, the in- 
variant mass is required to be within 3cr of the central 
value, where a is the fitted mass resolution. The invari- 



ant mass is then constrained to the nominal value |l(. 
Each resulting candidate passing the selection 

criteria is then combined with a positively charged track 
interpreted as a tt + (K + ) to form a S° candidate. For 
the f2 _ K + final state, the two tracks must be iden- 
tified as kaons. Particle identification is performed with 
dE/dx and Cherenkov angle measurements [Tlj . 

Additional selection criteria, described below, are used 
to improve the signal-to-background ratio. As a precau- 
tion against selection bias, these are optimized with sub- 
samples of the data: 20 fb" 1 and 40 fb -1 for the E~n + 
and S1~AT + final states, respectively. A minimum decay 
distance of 2.5 mm (1.5 mm) between the event primary 
vertex and the S - (fi - ) decay vertex in the plane per- 
pendicular to the beam direction is required. The dis- 
tance between the Q~ and A decay vertices is required 
to be at least 3 mm. In addition, the relative positioning 
of vertices is required to be causally connected: we re- 
ject candidates in which the S~ decays further from the 
primary vertex than its daughter A does, or where the 
displacement vector from the tt~ decay point to the A de- 
cay point is anti-parallel to the A momentum vector • 
The invariant mass distributions for the S° candidates 
in the full data set satisfying these criteria are shown in 
Fig. n (a) and (b) for S~7r + and £!~AT + combinations, 
with signal yields of approximately 8100 and 1000 events, 
respectively. 

Simulated events with the decaying into the two de- 
sired final states are generated for the processes e + e~ — > 
cc -> E° C X and e+e" -> T(4S") — > BB — > E° C X, where 
X represents the rest of the event. The pythia simu- 
lation package [l3|, tuned to the global BABAR data, is 
used for the cc fragmentation and for B decays to 
and GEANT4 0] is used to simulate the detector re- 
sponse. For cc production, samples of 90,000 events for 
the S~7r + final state and 60,000 for the £l~K + final state 
are generated. For BB production, samples of 255,000 
and 120,000 events are used, respectively. 

Additional generic Monte Carlo events are used to in- 
vestigate possible background contributions. The sample 
sizes are equivalent to 245 fb" 1 , 64 fb" 1 , and 33 fb" 1 for 
e + e~ — > BB, cc, and qq, respectively, where q = u,d, s. 
Excluding signal contributions, the mass distribution 
varies smoothly throughout the region near the S° mass. 

To measure the ratio of branching fractions, a further 
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FIG. 1: Invariant mass distributions for H° candidates in 
116.1 fb^ 1 of data, for (a) E~tt + , and (b) Q~ K+ . 



requirement that p* > 1.8 GcV/c is imposed on the 
S° candidates in order to suppress combinatoric back- 
ground and improve the signal purity. Additionally, the 
candidates are required to be within the region of high 
S° reconstruction efficiency —0.8 < cos 9* < 0.6, where 
9* is the polar angle of the S° candidate with respect 
to the collision axis in the center-of-mass frame. Af- 
ter these criteria, the signal yields for the E~n + and 
il~K + modes are approximately 3650 and 650, respec- 
tively. The efficiency is calculated from signal Monte 
Carlo events as a function of p* and cos#* for each of 
the decay modes. For each mode, a fifteen-parameter fit 
gives a smooth parameterization of the efficiency with 
small statistical uncertainty. The efficiency is then cor- 
rected by weighting each candidate by the inverse of its 
efficiency, and the efficiency-corrected mass spectrum is 
fitted with a double Gaussian with a common mean for 
signal plus a linear background function. Including effi- 
ciency loss due to the f2 _ and A branching fractions, we 
obtain 25889 ± 516 weighted events in the E~tt + mode 
and 7615 ±443 weighted events in the tt~K + mode. The 
X 2 fit probabilities are 65% and 5%, respectively. In each 
case, the wider Gaussian contributes approximately one 
quarter of the yield. 

We evaluate several sources of systematic uncertainty 
in the ratio of branching fractions: the fits to the mass 
spectra (3.4%), the efficiency parameterization (3.1%), 
particle identification (2.0%), finite Monte Carlo statis- 
tics (1.4%), multiple candidates in the same event (1.0%), 
charge asymmetries in detection efficiency (1.0%), the 



cos#* distribution (1.0%), and the f2~ branching fraction 
(1.0%). No baryon polarization is considered and any 
systematic uncertainty due to this is neglected. Adding 
all of the uncertainties in quadrature, we obtain: 



"0 _^ 77-„ 



= 0.294 ±0.018 ±0.016. 



After obtaining the ratio of branching fractions, we 
next measure the p* spectrum of the S° baryons in order 
to study the production mechanisms in both cc and BB 
events. The same selection criteria and data samples de- 
scribed above are used, except that no requirement on p* 
or cos 9* is made. Instead, the Sg candidates are divided 
into intervals of p* . The yield is then measured in each 
interval with two different methods: first with a fitting 
method, where the mass spectrum is fitted with a single 
Gaussian for signal plus a linear background function and 
the integral of the Gaussian is taken as the yield; second 
with a counting method, where the background is esti- 
mated from mass sidebands and the signal yield is then 
taken as the statistical excess above this background in a 
mass window around the peak. The use of two different 
methods serves as a cross-check. 

The efficiency in each p* interval is estimated with sig- 
nal Monte Carlo events from that p* range. For both 
methods, the simulated events are reconstructed and the 
yield is measured, then divided by the number of events 
generated to obtain the efficiency. Due to the different 
angular distributions, the efficiencies for Sg produced 
from cc (e c -) and from BB (s B g) differ slightly. In 
the region 1.2 < p* < 2.0 GeV/c where both produc- 
tion mechanisms are significant and the difference is ap- 
proximately 8% (relative), the efficiency is taken to be 
(e cE + e BB -)/2. The systematic uncertainty on the effi- 
ciency is then \e C c — £ BB \ 



|/V 12. The angular distribu- 



tions produced in pythia fragmentation are assumed to 
be correct when calculating the efficiency; the data are 
fully consistent with these distributions within available 
statistics. The efficiency-corrected yield in each p* inter- 
val is then calculated, including loss of efficiency due to 
the A and fl~ branching fractions. The spectra obtained 
with the two methods are in good agreement; we use the 
counting method for the quoted results since it is more 
stable for low statistics. 

A number of systematic uncertainties are considered, 
the most important of which are the uncertainties asso- 
ciated with the track-finding and particle identification 
efficiencies (5.8% and 3.5%, respectively). Uncertainties 
from the simulated S° mass resolution (1%), the mass 
resolutions of the intermediate hyperon states (0.5%), the 
p* resolution [(0(1%))], the effect of finite interval width 
[(0(2%))], multiple candidates (0%), non-linearity of the 
background [(0(1%))], the signal measurement method 
used (2%), the finite Monte Carlo statistics available 
[(0(3%))], and uncertainties in the A and f2~ branch- 
ing fractions (0.8%, 1.0%) are all considered individually; 
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FIG. 2: The p* spectrum measurements. In (a), the p* spec- 
trum of H° decaying via H~7r + is shown for the on-resonance 
data sample. In (b), the on-resonance and off-resonance data 
samples are shown together, with the off-resonance normal- 
ization scaled to account for the difference in integrated lu- 
minosity and cross-section. In each plot, the inner error bars 
give the statistical uncertainty and the outer error bars give 
the sum in quadrature of the statistical and systematic un- 
certainties. The vertical line at 2.15 GeV/c in (b) shows the 
kinematic cutoff for H° produced in B decays at BABAR.. Note 
that the vertical axes show events per unit p* , not events in 
each p* bin as given in Table [I] 



the notation 0(x%) indicates the typical value when the 
exact uncertainty varies among p* intervals. The total 
systematic uncertainty for each p* interval is obtained 
by adding the individual contributions in quadrature. In 
addition, a systematic correction of 1.0% is applied to ac- 
count for a known data-Monte Carlo discrepancy in the 
track-finding efficiency, and small corrections are applied 
to each interval to account for the broadening effect of 
the p* experimental resolution on the spectrum. The fi- 
nal p* spectrum for the on-resonance data set, obtained 
with the counting method in the 2 _ 7r + mode, is shown 
in Fig.^a). Table |U shows the corresponding values. 

A further check is performed by comparing the two 
decay modes. The Q,~K + yields are scaled by a factor 
of (1/0.294), the ratio of branching fractions previously 
presented in this Letter. Because the il~K + signal has 
fewer events, wider p* intervals are used. The spectra 
of the two modes show good agreement in both shape 
and normalization and have a \ 2 probability of 80% for 
consistency. This serves as a cross-check both of the p* 
spectrum measurement and of the ratio of branching frac- 
tions. 

The double-peak structure seen in the p* spectrum is 



TABLE I: Efficiency-corrected yield and cross-section prod- 
uct including B production a(e + e~ — > H°A") x B{E.° c — + 
H~7T + ), for the on-resonance data. 



p* range 


Corrected 


Cross-section 


(GeV/c) 


yield 


product (fb) 


0.0- 


-0.2 


1046 ± 201 ± 128 


10 ± 2 ± 1 


0.2- 


-0.4 


5889 ± 446 ± 483 


56 ± 4 ± 5 


0.4- 


-0.6 


10681 ± 631 ± 801 


101 ± 6 ± 8 


0.6- 


-0.8 


10709 ± 660 ± 817 


102 ± 6 ± 8 


0.8- 


-1.0 


8811 ± 647 ± 679 


84 ± 6 ± 7 


1.0- 


-1.2 


6834 ± 573 ± 530 


65 ± 5 ± 5 


1.2- 


-1.4 


2954 ± 501 ± 252 


28 ± 5 ± 2 


1.4- 


-1.6 


2429 ± 470 ± 212 


23 ± 4 ± 2 


1.6- 


-1.8 


2252 ± 424 ± 202 


21 ± 4 ± 2 


1.8- 


-2.0 


2159 ± 350 ± 217 


20 ± 3 ± 2 


2.0- 


-2.2 


2375 ± 347 ± 205 


23 ± 3 ± 2 


2.2- 


-2.4 


2743 ± 340 ± 227 


26 ± 3 ± 2 


2.4- 


-2.6 


3537 ± 315 ± 285 


34 ± 3 ± 3 


2.6- 


-2.8 


3920 ± 282 ± 306 


37 ± 3 ± 3 


2.8- 


-3.0 


4595 ± 294 ± 359 


44 ± 3 ± 3 


3.0- 


-3.2 


4873 ± 263 ± 401 


46 ± 2 ± 4 


3.2- 


-3.4 


4449 4- 944 -1- ^48 

4:4:4:^ ZIZ jl^i^ ZH 04:0 


42 ± 2 ± 3 


3.4- 


-3.6 


4084 ± 223 ± 355 


39 ± 2 ± 3 


3.6- 


-3.8 


2282 ± 171 ± 189 


22 ± 2 ± 2 


3.8- 


-4.0 


2095 ± 155 ± 166 


20 ± 1 ± 2 


4.0- 


-4.2 


1168 ± 123 ± 177 


11 ± 1 ± 2 


4.2- 


-4.4 


233 ± 53 ± 32 


2.2 ± 0.5 ± 0.3 


4.4- 


-4.6 


88 ± 37 ± 21 


0.8 ±0.3 ±0.2 


4.6- 


-4.8 


5 ±13 ±7 


0.0 ±0.1 ±0.1 


4.8- 


-5.0 


24 ± 17 ± 16 


0.2 ± 0.2 ± 0.1 



TABLE II: Cross-section product including B production 
a{e+e- -» E° C X) x B(S° C -> S~ir + ), for the on- and off- 
resonance data. The off-resonance cross-sections are scaled 
to a center-of-mass energy of 10.58 GeV. 



p range 
(GeV/c) 



Cross-section 
On-resonance 



product (fb) 
Off-resonance 



0.00- 


-0.45 


88 ± 5 ± 7 


10 ± 12 ± 1 


0.45- 


-0.90 


218 ±9± 17 


16 ±21 ±2 


0.90- 


-1.35 


128 ±8± 10 


-7 ±20 ±2 


1.35- 


-1.80 


51 ±6 ±4 


16 ±18 ±2 


1.80- 


-2.15 


37 ± 4 ± 3 


23 ±13 ±2 


2.15- 


-2.70 


83 ± 5 ± 6 


91 ±16 ±7 


2.70- 


-3.30 


133 ±4± 10 


168 ± 15 ± 13 


3.30- 


-4.00 


99 ± 3 ± 8 


89 ± 10 ± 7 


4.00- 


-4.70 


14 ± 1 ± 1 


17±4±2 
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due to two production mechanisms: the peak at lower 
p* is due to S° production in B meson decays and the 
peak at higher p* is due to E° c production from the cc 
continuum. This is evident from Fig.^b), where the p* 
spectra for the on-resonance and off-resonance data are 
shown separately (with the off-resonance spectrum scaled 
to the on-resonance integrated luminosity and corrected 
for the change in cc cross-section) . Table [H] shows the 
corresponding values. The cc peak is present in both 
samples, but the BB peak is only present in the on- 
resonance sample. Assuming baryon number conserva- 
tion, the kinematic limit for E° c produced in the decays 
of B mesons at BABAR is p* = 2.135 GeV/c. We com- 
pare the on-resonance and scaled off-resonance samples 
for p* < 2.15 GeV/c to obtain the yield of E° c produced 
in B decays. This is scaled by the number of B mesons in 
the data sample (introducing a further 1.1% systematic 
uncertainty) to obtain: 

B(B E° C X) x B(E° C S-7T+) 

= (2.11 ±0.19 ±0.25) x 10~ 4 . 

The yield of produced in cc events at an energy of 
10.58 GeV is calculated from the scaled off-resonance 
data set (for p* < 2.15 GeV/c) and the on-resonance 
data set (for p* > 2.15 GeV/c). The yield is then divided 
by the integrated luminosity (introducing a further 1.5% 
systematic uncertainty) to obtain the cross-section from 
the continuum: 

a(e+e- -> E° C X) x B(E° C -> E'n+j 
= (388 ±39 ±41) fb, 

where both and E° c are included in the cross-section. 
The effect of initial state radiation is not isolated. 

In summary, we have studied the S° charmed baryon at 
BABAR through its decays to the fl~K + and E~ir + final 
states using 116.1 fb -1 of data. The ratio of branch- 
ing fractions of these decay modes was measured to be 
0.294 ± 0.018 ± 0.016. This represents a substantial im- 
provement on the previous measurement |3( and is con- 
sistent with a quark model prediction We have also 
measured the p* spectrum for E° c produced at the T(4S*) 
resonance. The high rate of E° c production at low p* 
in B decays (below 1.2 GeV/c) is particularly intrigu- 
ing, implying that the invariant mass of the recoiling an- 
tibaryon system is typically above 2.0 GeV/c 2 . This can 
be explained naturally by a substantial rate of charmed 
baryon pair production through the b — > ccs weak decay 
process [a, [jj, |8j, |9( which was observed indirectly in a pre- 
vious BABAR analysis 0|. In this Letter we measured 



the branching fraction product B(B — > x B(z? c — > 

E-tt + ) to be (2.11 ± 0.19 ± 0.25) x 10~ 4 ; the preci- 
sion i s sig nificantly improved over the previous measure- 
ment |Uj . We have also measured the cross-section prod- 
uct cr(e + e~ — * S°A) x B{E a c — ► E~ir + ) from the contin- 
uum to be (388 ± 39 ± 41) fb. 
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